Nonequilibrium spin-dependent phenomena in mesoscopic superconductor-normal 

metal tunnel structures 
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We analyze the broad range of spin-dependent nonequilibrium transport properties of hybrid 
systems composed of a normal region tunnel coupled to two superconductors with exchange fields 
induced by the proximity to thin ferromagnetic layers and highlight its functionalities. By calculat- 
ing the quasiparticle distribution functions in the normal region we find that they are spin-dependent 
and strongly sensitive to the relative angle between exchange fields in the two superconductors. The 
impact of inelastic collisions on their properties is addressed. As a result, the electric current flow- 
ing through the system is found to be strongly dependent on the relative angle between exchange 
fields, giving rise to a huge value of magnetoresistance. Moreover, the current presents a complete 
spin-polarization in a wide range of bias voltages, even in the quasiequilibrium case. In the nonequi- 
librium limit we parametrize the distributions with an "effective" temperature, which turns out to 
be strongly spin-dependent, though quite sensitive to inelastic collisions. By tunnel coupling the 
normal region to an additional superconducting electrode we show that it is possible to implement 
a spin-polarized current source of both spin species, depending on the bias voltages applied. 

PACS numbers: 72.25.-b,85.75.-d,74.50.+r,05.70.Ln 



I. INTRODUCTION 

Although the interest in nonequilibrium superconduc- 
tivity dates back to the seventies^, nonequilibrium trans- 
port phenomena in hybrid superconducting structures 
are currently under the spotlight. One of the key ex- 
periments that renewed this interest was probably the 
control of the supercurrent flowing through a Josephson 
junction, and even the reversal of its sign, accessible by 
altering the quasiparticle population in the weak link (see 
Ref. [U and references therein). Out-of-equilibrium elec- 
tron population can be realized in mesoscopic conduc- 
tors subject to a bias voltage in which electrons cannot 
exchange energy either with one another or with lattice 
phonons, so that their energy distribution is not Fermi 
like 3 -. Quasiequilibrium is reached if electrons can ther- 
malize, while still decoupled from the phonons, so that 
they can reach a temperature which is different from the 
one relative to the phonon bath. In ballistic Josephson 
junctions supercurrent control occurs by inducing a non- 
equilibrium population of Andreev levels either by in- 
jecting a current through an additional normal terminal 
connected to the weak link^ or by applying an elec- 
tromagnetic radiation on the weak link^i. The diffu- 
sive long-junction limit was considere d 8 ' 9 ' 1 ' 1 1 and ex- 
perimentally realized too 12 ' 13 . The control of supercur- 
rent by cooling electrons in the weak link was proposed 

m Refs. EmiiKn and experimentally realized^. It is 

worthwhile stressing that electron temperature can be 
lowered below the phonon temperature, thus realizing 
electron microrefrigeratio n 19 ' 20 , by explo iting the super- 
conducting energy gap (see Refs. |2l|22l 

and references 

therein). 

Spin-dependent properties in out-of-equilibrium hy- 
brid systems were investigated in a limited number of ar- 



ticles. In Refs. ferromagnet-superconductor- 
ferromagnet (FSF) double tunnel junctions were consid- 
ered in order to study the spin imbalance induced in S 
by non-equilibrium. In the anti-ferromagnetic alignment 
of the magnetizations of the F layers a strong suppres- 
sion of superconductivity was found, leading to a large 
magnetoresistive effect. In Josephson junctions the effect 
of spin injection^! and presence of weak ferromagnets 28 
was considered, while the effect of Andreev reflection on 
spin accumulation in a ferromagnetic wire was reported 
in Ref. [29|. In Ref. [3(3 the possibility of manipulating mag- 
netism through the interplay of superconductivity and 
nonequilibrium transport was investigated. Recently we 
have proposed 3 ^ a hybrid ferromagnet-superconductor 
(FS) spin valve whose operation is based on the inter- 
play between out-of-equilibrium quasiparticle dynamics 
and proximity-induced exchange coupling in supercon- 
ductors. Huge tunnel magnetoresistance values as high 
as several 10 6 % has been predicted, leading to a fully- 
tunable structure which shows high potential for appli- 
cation in spintronics. In this paper we comprehensively 
investigate the physics and functionality of the setup an- 
alyzed in Ref. [3l|, extending our study to the presence of 
finite electron-electron interaction and to the quasiequi- 
librium limit, as well as to the presence of nonidealities 
in the superconductors. In this setup a spin-dependent 
"effective" temperature for the electrons in the N region 
emerges, thus leading to possible new spin-dependent 
thermoelectric effects. 

The paper is organized as follows: in Sec. [TT] we de- 
scribe the system under investigation and in Sec. IIIII we 
derive the quasiparticle distribution functions in different 
regimes. In particular, we consider the nonequilibrium 
limit in Sec. IHI Al we include the effect of inelastic colli- 
sions in Sec. IIII B[ and we describe the quasiequilibrium 
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regime in Sec. MI CI In Sec. IIVI we discuss the behavior 
of the electric current, focusing on the magnetoresistive 
effects and on the spin-filtering properties of the system 
in Sec. [V] Section I VII is devoted to the characterization 
of the nonequilibrium distribution through an "effective " 
temperature, and to the exploitation of the system as a 
source of spin-polarized current through the introduction 
of an additional superconducting electrode. Finally, we 
draw our conclusions in Sec. IVIIl 



II. SETUP 

We consider a device consisting of two identical FS bi- 
layers (FSi 2) symmetrically connected to a mesoscopic 
normal metal region (N) of length through tunnel 
contacts (I) of resistance i?t- The concentration of im- 
purities is such that quasiparticle transport is diffusive. 
The resulting system, a FS-I-N-I-SF heterostructure, is 
shown in Fig. [T] in two distinct experimental implemen- 
tations. Figure [IJa) shows a spin valve-like structure, 
which consists of a sequence of stacked metallic layers, 
while Fig. [1Tb) displays a planar system. Although the 
two implementations are equivalent on theoretical foot- 
ing, the planar configuration allows the measurement of 
local properties (e.g., the quasiparticle distribution func- 
tions as well as the local temperature) by connecting the 
N region to additional metallic probes. This will be ad- 
dressed in Sec. IVIl For the sake of simplicity we assume 
a symmetric system (a resistance asymmetry would not 
change the overall physical picture), t F (ts) labels the F 
(S) layer thickness and a bias voltage V is applied across 
the structure. The exchange field in the left ferromagnet 
(hi) is aligned along the z axis for the setup in Fig. [JJa) 
or along the y axis for the setup in Fig. UJb), while that 
in the right F layer (h 2 ) is misaligned by an angle <j> [see 
Fig. DTa'),(b')]. For simplicity we set \hi\ = \h 2 \ = h. 
In real structures h 2 can be rotated by applying an in- 
plane magnetic field as low as some mT. Moreover we 
assume that (i) the FS interface is transparent and (ii) 
Rt is much larger than both the resistance of the N layer 
(.Rn) and the FS contact resistance. The first condition 
ensures that the superconductor is strongly affected by 
the proximity of the F layer—, while the second ensures 
that all the voltage drop occurs at the tunnel barriers 
(so that any spatial variation of the chemical potential 
within the N region can be neglected), and that each FS 
bilayer is in local equilibrium. 

The electronic properties of a FS bilayer can be 
analyzed within the quasiclassical Green's function 
formalism^. We are interested in the situation in which 
the influence of the F layer on the superconductor be- 
comes nonlocal. This occurs in the limit tg < £g = 
y / hD/2nkBT c and t F < £p = y/hD/h, where £3 and £f 
are the superconducting coherence length and the length 
of condensate penetration into the ferromagnet, respec- 
tively. D denotes the diffusion coefficient, T c is the super- 
conducting critical temperature and fee is the Boltzmann 




FIG. 1: (color online) Sketch of two possible implementa- 
tions of the FS-I-N-I-SF structure analyzed in this work, (a) 
Spin-valve setup consisting of a sequence of stacked metal- 
lic layers, (b) A planar structure. Ferromagnetic layers (F) 
induce in each superconductor, through the proximity effect, 
an exchange field (/ii,a) whose relative orientation can be 
controlled by an externally applied magnetic field. The F 
exchange fields are confined (a') to the y — z plane for the 
setup shown in (a), and (b') to the x — y plane for the setup 
shown in (b), and are misaligned by an angle <f>. A voltage 
bias V, applied across the structure, allows to control the en- 
ergy distributions in the N region. The structure is assumed 
quasi-one-dimensional. 



constant. In this situation, the ferromagnet induces in 
S a homogeneous effective exchange field (analogous to 
the one present in magnetic superconductors^) through 
proximity effect and modifies the superconducting gap 
(A). The effective values of the exchange field (h*) and 
gap (A*) are given by2£: 

A*/A = v s t s (v s t s + v F t F )- 1 

h*/h = v F t F (vsts + VFt F )- 1 , { ' 

where v§ (v F ) is the normal-state density of states (DOS) 
in S (F). fn particular, if v F = and for t F <C ts, it 
follows that 



h*/h ~ tp/ts < f, y ' 

i.e., h* turns out to be much smaller than in an iso- 
lated F layer. As a matter of fact, h* can take values of 
the order of magnitude of A*. These conditions can be 
achieved quite easily in a realistic structure. We assume 
that the only effect of h* on the quasiparticles is to lead 
to a spin-dependent superconducting DOS, i.e., we ne- 
glect any influence of the induced magnetic moment on 
the orbital motion of electrons. Furthermore, we assume 
negligible spin-orbit interaction^. The superconductor 
DOS (A/f ) thus will be BCS-like, but shifted by the ef- 
fective exchange energy (equivalent to that of a Zeeman- 
split superconductor in a magnetic field^). By choosing 
the spin quantization axis along the direction of the ex- 
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FIG. 2: (color online) Spin-dependent quasiparticle distribu- 
tion functions fa{e) in the full nonequilibrium limit vs en- 
ergy e for several bias voltages at <j> = 0, Tb at h = 0.1T C , and 
ft* =0.2 A*, (a) / + ( e ); (b) /_(e). 



change field we have 



Re 



e + ah* + iT 



y/(e + ah* + IT) 2 



A* 



. (3) 



where e is the energy measured from the conden- 
sate chemical potential, a = ±1 refers to spin par- 
allel (antiparallel) to the direction of h%, and T is a 
smearing parameter—. The latter allows quasiparticle 
states within the gap due to inelastic scattering in the 
superconductor 3 -, or inverse proximity effect from the 
nearby metallic layers. Typical values for T lie in the 
range T ~ 1 x 10~ 5 A . . . 1 x 10~ 3 A for Al as a thin-film 
superconducting electrode^. In the following calcula- 
tions we set r = 10 _4 A*, unless differently stated. 

In order to describe our system we make use of the 
tunneling Hamiltonian approach, and neglect proximity 
effects at NIS interfaces. 



III. QUASIPARTICLE DISTRIBUTIONS 



Negligible inelastic scattering: 
nonequilibrium limit 



full 



At finite bias V and in the limit of negligible inelas- 
tic scattering, quasiparticles in the N layer will be out of 
equilibrium and thus, in general, not distributed accord- 
ing to the Fermi function. The steady-state nonequilib- 
rium distribution functions can be calculated by equat- 
ing, at each energy value, the tunneling rate of quasi- 
particle entering the N region from the insulating layer 
on the left-hand-side to the tunneling rate of those ex- 
iting through the right-hand-side barrier—. In the gen- 
eral case of non-collinear exchange fields, the spin eigen- 
states relative to S2 (| t) an d I I)) can be obtained by 



rotating the spin eigenstates relative to Si (|+) and |— )) 
by the angle cf> (representing the misalignment between 
hi and h^)- As a consequence, spin up {a = +1 with 
eigenstate |+)) quasiparticles exiting the N layer through 
the right-hand-side barrier will now consist of two con- 
tributions. One describes tunneling into spin up (with 
eigenstate | f)) quasiparticles, proportional to cos 2 [0/2], 
and the other describing tunneling into spin down (with 
eigenstate | [)) quasiparticles, proportional to sin 2 [0/2]. 
As a result, the nonequilibrium distribution function in 
the N layer is spin-dependent and can be written as 



fa(e,V,h* 



Na 1 + a(4>)Af(T 2 + K<f)N% 

(4) 

where a{<f) = cos 2 [0/2], b(<j>) = sin 2 [0/2], T Sl( - s ^ = 
fo(e±eV/2), A/^ = ^(e + eV/2), A^ 2 = ^(e-eV/2), 
fo(s) is the Fermi function at bath temperature Tbath, 
and e is the electron charge. 

Figures EJa) and (b) show the nonequilibrium dis- 
tributions functions [calculated from Eq. (J4j)] for spin 
up and spin down quasiparticles, respectively, vs en- 
ergy e for the parallel configuration (i.e., <j> — 0) at 
T bath = 0.1T C , h* = 0.2A* and different values of V (we 
assume the superconducting gap to follow the BCS rela- 
tion A* = 1.764fceT c ). Figure[5]shows that, by increasing 
the bias voltage V , spin up and spin down distributions 
are shifted in opposite directions on the energy axis, sim- 
ilarly to what is expected in the presence of an effective 
spin-dependent chemical potential (fJ^). I n particular 
f+(e) is shifted toward negative energies, while /-(e) to- 
ward positive energies. Moreover, for eV > A*, the spin- 
dependent chemical potential saturates at /i° ff = —ah*. 
As shown in Ref . [3(J this effect can be used to electrostat- 
ically manipulate the magnetic properties of the N region. 
The role of a finite T (i.e., the presence of quasiparticle 
states within the gap) can be appreciated in Fig. O By 
increasing eV from to A* the distributions broaden, 
and reflect the "heating" of the N region, as discussed 
in Refs. [HH This effect is absent for T = 0. By fur- 
ther increasing the bias voltage the distribution functions 
sharpen due to "cooling" provided by the superconduct- 
ing energy gap^ 2 -. 

Analogously, in Figs.[3{a) and (b) we plot the nonequi- 
librium distribution functions for spin up and spin down 
quasiparticles, respectively, for the antiparallel configu- 
ration (i.e., 4> = 7r). Distribution functions are shown vs 
energy e for different values of V, and were calculated 
for the same parameters as in Fig. O In this case, up 
and down distributions remain centered around e = 
upon biasing (equivalently, their effective chemical po- 
tential is always /i° ff = 0), but at a given bias voltage 
the features of the distributions are more pronounced for 
the spin down case. As we shall see in Sec. VI, up and 
down distributions are characterized by different effective 
electronic temperatures (T° ff ). In general, for any angle 
4> differing from or 7r the spin-dependent distribution 
functions f a (e) will be characterized by both an effective 
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FIG. 3: (color online) Spin-dependent quasiparticle distribu- 
tion functions fa{e) in the full nonequilibrium limit vs en- 
ergy e for several bias voltages at 4> — ir, Tbath = 0.1T C , and 
h* = 0.2A*. (a) /+(£); (b) 
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FIG. 4: (color online) Spin-dependent quasiparticle distribu- 
tion functions fcr(e) vs energy e calculated for several /C co n 
values at 4> = 0, eV = A*, T ba th = 0.1T C , and h* = 0.2A*. 
(a) / + ( e ); (b) /_( e ). 



chemical potential and an effective electronic tempera- 
ture. 



B. Intermediate inelastic scattering 



tions functions f a , denned by 



In the presence of scattering the approach of Sec. IIII Al 
cannot be used and one has to resort to the kinetic equa- 
tion theory. Electrons in metals experience both elas- 
tic and inelastic collisions. The latter drive the sys- 
tem to equilibrium and can be expected to hinder the 
manifestation of the phenomena discussed in the previ- 
ous section. At low temperatures (typically below 1 K) 
electron-electron scattering 39 and scattering with mag- 
netic impuritie a 40 i 41 are the dominant sources of inelastic 
collisions ^ 41 ' 42 . Since Rt is in general large compared to 
the wire resistance [R^ = tt;/(Afpe 2 DA)], where Afp is 
the N-region DOS at the Fermi energy and A the wire 
cross-section), we can assume that f a does not depend 
on the position in the wire 1 ^. 

In the following we shall analyze the role of inelastic 
electron-electron relaxation on the quasiparticle distri- 
bution. The effect of electron-electron scattering due to 
Coulomb interaction on the spin-dependent distributions 
can be accounted for by solving a pair of coupled station- 
ary kinetic equations: 



D 



dx 2 
d 2 /-(g) 
dx 2 



(5) 



J coll( £ )' 



together with the Kuprianov-Lukichev boundary condi- 
tions at the NIS interfaces 4 ^. In Eqs. (© 2^ oll (e) is the 
net collision rate at energy e, functional of the distribu- 



(6) 



where 



fc(w) 



dE {f+(E + u) [1 - /+ (E)} + f-(E + u) [1 - /_ (E)]} 

(7) 



and 



J dE{f + (E)[l-f + (E- 



^[l-« e - W )] 



f_(E)[l-f_(E + u)}}. 

(8) 



In Eqs. (0 and |8]), k(co) = k cc uj 3 / 2 according to the 
theory of screened Coulomb interaction 4 ^ for a quasi-one 



^2Dh 3 / 2 AT^A) 



-1 45,46 
15 



dimensional wire, where k cc = (- 
By rewriting Eqs. ([5]) in dimensionlcss unit 
the strength of the electron-electron interaction can 
be expressed as /C co u = (R t / R^it^n^/ 'D)y/K = 
(t N /V^)(R t /R K )y/A/hD, where R K = h/2e 2 . We note 
that the strength of the electron-electron interaction 
turns out to be proportional to the length of the wire 
as well as to the tunnel barrier resistance. 

We solved Eqs. © with h* = 0.2 A*, eV = A* and 
T = 0.1T C for several /C co ii values 4 ^. The effect of 
electron-electron scattering on the quasiparticle distri- 
bution functions is displayed in Figs. [4] and [5] for (f> = 
and (f> = it, respectively. For the <fi — case, by in- 
creasing Coulomb interactions the quasiparticle distribu- 
tions are forced toward thermal ones still characterized 
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FIG. 5: (color online) Spin-dependent quasiparticle distribu- 
tion functions / CT (e) vs energy e calculated for several /C co n 
values at <f) = n, eV = A*, T ba th = 0.1T C , and h* = 0.2 A*, 
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by different chemical potentials for both spin species^. 
In the antiparallel configuration (see Fig. O the ef- 
fect of inelastic relaxation is similar, but now the spin- 
dependent distribution function will coincide for suffi- 
ciently large /C C oii values. It is easy to recognize that, in 
both thermal Fermi-like distribution is reached 

for tC co \\ of the order of 10. Assuming parameters for 
a realistic Al/Al 2 3 /Ag SINIS microstructur o 21 ' 36 (with 
A ~ 200 ^eV, D = 0.02 m^" 1 and R t = 1 kfi) JC coU = 10 
corresponds to a rather long N region, in « 47 /mi. 



Strong inelastic scattering: quasiequilibrium 
limit 



This is the regime characterized by the fact that the 
electron-electron interaction is so strong that quasipar- 
ticles can reach an equilibrium (Fermi-like) distribution, 
while the electron-phonon coupling is negligible 2 ^. Such 
distributions are characterized by quasiequilibrium chem- 
ical potential and temperature. Since electron-electron 
interaction occurs between quasiparticles irrespective of 
their spin (in the absence of spin- mixing mechanisms), 
the quasiequilibrium temperature (T qe ) will be indepen- 
dent of spin, and different from the temperature of the 
phonon bath Tbath- On the contrary, since electron- 
electron interaction redistributes the energy among elec- 
trons of a given spin species, in the absence of spin- 
mixing mechanisms, the quasiequilibrium chemical po- 
tential (fi^ c ) will depend on spin. This is a consequence 
of the fact that the number of electrons of a given spin 
must be conserved. Furthermore, both quasiequilibrium 
chemical potential and quasiequilibrium temperature will 
depend on 0, therefore they will be different for parallel 
and antiparallel configurations. 



In the absence of spin-flip mechanisms, the quasiequi- 
librium distribution functions can be calculated by im- 
posing the conservation of particle currents, indepen- 
dently for the two spin species, together with a balance 
equation for the heat currents. In particular, in the for- 
mer case we require that 



%(v,h*,<t>) = i?(V,h*,4>), 



(9) 



where 
and 



1 

7R t 



de^(e)[^(e)-U(e)} (10) 



J*(V, h*A)= ^-Jde (e)o(^) + M s Ue)b(0) 



11) 



are the electric currents flowing through the left/right 
(L/R) NIS interface. Note that, in contrast to the full 
nonequilibrium regime where the the tunneling rates are 
set to be equal at each energy, here the conservation in- 
volves the total currents, since the electron-electron inter- 
action mixes the energy of the electrons. In the absence 
of electron-phonon coupling, the only contribution to the 
heat flux is the heat current flowing off the N region 
through each NIS interface. The latter is given by 



J^V,h*,<f>) = -^-J dee (e) [f a (e) - ^ (e) 



for the left NIS contact, and by 



(12) 



jf{V,h* 



^-Jdee [N^ (e)a(<j>) + Af s Ue)b(cb) 

x[f a (s)-^(e)] (13) 



for the right contact. The balance equation for the heat 
flux thus simply reads 



E 



[jZ(v,h* t 4>) + j?(y t h* 



o. 



(14) 



By assuming that f a = fo(e — /Lt qe , T qe ) and solving © 
and (|14[) . the temperature and chemical potentials can be 
easily determined. It turns out that, while in the antipar- 
allel alignment spin up and down distributions are equal, 
in the parallel one the two spin components have equal 
effective electronic temperature (though different from 
the antiparallel alignment), but opposite effective chem- 
ical potential (see Figs. 0] and [5] for large K co \\ values). 
Although the quasiequilibrium regime might seem an un- 
realistic limit, it actually describes the case of strong 
electron-electron interaction quite well. Indeed, accord- 
ing to our calculations (Sec. III.B), quasiequilibrium dis- 
tributions are already reached for an electron-electron 
collision strength JC co \\ ~ 10. In the following Sections 
we shall investigate the impact of quasiequilibrium on 
spin-dependent transport properties. 
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IV. ELECTRIC CURRENT 

The transport properties of the FS-I-N-SF structure 
are determined by the spin-dependent distribution func- 
tions f a . We note that although a Josephson current 
can flow through the system, its theoretical description 
is beyond the scope of the present paper. As a mat- 
ter of fact, we shall be only concerned with quasiparticle 
transport. Furthermore, although similar results for tun- 
nel magnetoresistance and current polarization could be 
obtained in a FS-I-SF structure (i.e., without the N in- 
terlayer) and not relying on noncquilibrium, the present 
system possesses a crucial advantage. In fact a FS-I-SF 
structure implies an additional undesired Josephson cur- 
rent, which can be fairly large as compared to the quasi- 
particle current (around one order of magnitude larger 
than the quasiparticle current relevant for high tunnel 
magnetoresistance and current polarization, see for ex- 
ample Ref. l33l ). Such supercurrent could be suppressed, 
for instance, by the application of an additional in-plane 
magnetic field. This field, however, would largely exceed 
that required to control the orientation of h*. By con- 
trast, in the FS-I-N-I-SF system the supercurrent can 
be kept extremely small up to a large extent, depending 
mainly on t^, on the tunnel barriers transmissivity, and 
on the N-interlayer material parameters. A simple esti- 
mate for the Josephson coupling in our structure reveals 
that the supercurrent can be from one to several orders 
of magnitude smaller than the quasiparticle current (see, 
for example, Ref. 

The quasiparticle current I (e.g., evaluated at the left 
interface) is given by 



(15) 



Figure [f^a) displays the electric current in full noncqui- 
librium vs bias voltage V calculated for several angles 4> 
at h* = 0.2A* and T = 0.1T C . A sizable current starts 
to flow only when the voltage V is such that the DOS 
is finite for both superconductors in some range of ener- 
gies. For <f> = 0, the current rises sharply at \eV\ — 2A*, 
similarly to the quasiparticle current of a SIS junction 
(also in the presence of an in-plane magnetic field^). In 
this case, in fact, the DOS of a given spin is shifted by 
the Zeeman energy in the same direction for both super- 
conductors. In contrast, for <p = n current sets off at 
\eV\ = 2(A* - h*). 

Figure [6jb) shows the nonequilibrium differential con- 
ductance 



G(V, h* 



dI(V, h* 



dV 



(16) 



calculated for the same values as in Fig. [U(a). Ad- 
ditional features are present at \eV\ = 2h* which are 
strongly temperature-dependent, and vanish in the limit 
T — > (the zero-bias conductance peak for <f> 7^ tt re- 
sembles that typical of a SIS junction composed of iden- 
tical superconductors^). These are a consequence of 
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FIG. 6: (color online) (a) Nonequilibrium current vs bias volt- 
age V for several angles cj> at Tbath = O.lTc and h* — 0.2 A*, 
(b) Nonequilibrium differential conductance G vs V calcu- 
lated for the same values as in (a). 



the overlapping of the superconducting DOSs where only 
thermally-activated quasiparticles exist at finite temper- 
ature. 

All this simply reflects how the spin-dependent DOS 
in each superconductor contributes to the total quasi- 
particle current at different V. This can be easily visu- 
alized by inspecting Fig. which shows idealized finite- 
temperature exchange field-split superconducting DOS 
for parallel spin species, at different bias voltage V and 
for the case (f> = it. In this case the DOS of Si is shifted 
in the opposite direction with respect to that of S2 [see 
Fig. 0Ja) for V = 0]. Then, by biasing the structure, 
the required voltage for a current to flow is smaller with 
respect to the = case, i.e., eV = 2 (A* — h*) [see 
Fig. |7{b)]. In the same way, for negative voltages, the 
current sets off at eV — — 2(A* + h*), as shown in Fig. 
[7jd). It is also clear that antiparallel spin species will 
give rise to features at the opposite bias voltage, there- 
fore explaining the origin of additional feature appearing 
at eV = 2(A* + h*). For intermediate values of </>, fea- 
tures are present at \eV\ = 2(A*±h*) and at \eV\ = 2A*, 
since contributions from both = and <f> = 7r configu- 
rations are present. Of particular relevance is the voltage 
interval 2(A* — h*) < \eV\ < 2A*. By increasing <j) from 
to 7T, the current is enhanced from a vanishingly small 
value up to a finite value leading to a spin-valve effect. 
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FIG. 7: (color online) Idealized finite-temperature exchange 
field-split density of states 7V Sl ' 2 of Si, 2 for parallel spin 
species, at (j> = 7T and different bias voltage V. In partic- 
ular, (b), (c) and (d) show how features in the tunneling 
current originate at eV = 2(A* — h*), eV = —2h*, and 
eV = — 2(A* + h*), respectively. Antiparallel spin species 
gives rise to features at opposite voltages. Green-dashed lines 
represents the superconducting DOS in the absence of the 
exchange field. 



It is noteworthy to mention that the nonequilibrium 
condition is essential for the observation of the spin- valve 
effect. At equilibrium the distribution functions in the N 
layer would be thermal and spin- independent. 



V. MAGNETORESISTANCE 

The spin- valve properties of the FS-I-N-I-SF setup can 
be evaluated quantitatively by analyzing the tunnel mag- 
netoresistance ratio (TMR), defined as 



TMR(V, h* 



G(V,h*,(/>)-G(V,h*,0) 
G(V,h*,0) ' 



(17) 



Figure Ufa) displays the absolute value of the nonequilib- 
rium TMR vs bias voltage V calculated for several angles 
0atT = O.lT c and/i* = 0.2A*. For2(A*-/i*) < \eV\ < 
2A* the TMR increases monotonically by increasing <f> 
and is maximized at <j) = tt where it reaches huge values 
exceeding 10 6 %. We note that in the limit T = and 
r = 0, I TMR I diverges, realizing an ideal full spin- valve 
effect. The nonequilibrium TMR behavior for several ex- 
change field values is shown in Fig. Efb), at T = 0.1T C 
and <j) = 7r. By decreasing h* , the maximum TMR value 
reduces, and so does the voltage interval of larger mag- 
netoresistance. Larger h* values are thus preferable in 




FIG. 8: (color online) (a) Nonequilibrium tunnel magnetore- 
sistance ratio |TMR| vs V calculated for several angles <f) at 
Tbath = O.lTc and h* = 0.2A*. (b) |TMR| vs V for different 
h* values at T = 0.1T C and d> = tt. 



order to extend the voltage window for optimized opera- 
tion and to maximize the TMR. 

The spin-filtering properties of this system can be 
quantified by inspecting the current polarization (Pi), 
denned as 



Pi(V,h* 



i\(v,h\4>)-iHv,h* 



(18) 



The calculated nonequilibrium Pj vs V is displayed in 
Fig. |5fa) for several (j) values, at T = 0.1T C and h* = 
0.2 A*. Upon increasing <f), two intervals of 100% spin- 
polarized current develop for 2 (A* — h*) < \eV\ < 2A*, 
extending to wider regions [2(A* — h*) < \eV\ < 2(A* + 
h*)} as (f> approaches tt. For 4> = 0, Pj vanishes like in SIS 
junctions with an in-plane magnetic field 35 . Depending 
on bias, fully spin-polarized currents of both parallel and 
antiparallel spin species can be obtained. The structure 
can thus be also operated as a controllable spin-filter by 
changing the orientation of h 2 as well as by varying V. 
Figure [Hb) shows Pi vs V for several h* at T = 0.1T C 
and <f> = tt. The net effect of increasing h* is to widen 
the regions of 100% spin-polarized current. 

It is important to discuss the effect of the smearing 
parameter T (which controls the presence of quasiparti- 
cle states within the superconducting gap) on the mag- 
netoresistance and current polarization. As shown in 
Fig. [TUT a). by increasing T, the TMR value decreases 
mostly in the region 2(A* - h*) < \eV\ < 2A*, while 
for other values of V almost no changes are found apart 
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FIG. 9: (color online) (a) Nonequilibrium current polariza- 
tion Pi vs V calculated for several angles <f> at Tbath = 0.1T C 
and h* = 0.2A*. (b) Pi vs V for different h* values at 
T bath = 0.1T c and = tt. 



from some smoothing of sharp features. In particular, 
the normal character of transport is strengthen by in- 
creasing r which causes a suppression of the large TMR 
value. The latter indeed is a consequence of the presence 
of the superconducting gap. On the contrary, the impact 
of r on Pi, plotted in Fig. [POT b) as a function of the 
voltage V, is much weaker: the polarization in the range 
2(A* - h*) < \eV\ < 2(A* + h*) is almost insensitive 
to r, being slightly reduced only for T values as large as 
1(T 2 A*. 

TMR values are expected to be marginally affected 
by the presence of electron-electron relaxation in the 
N layer. Indeed, as discussed in Sees. IIIIBI and IIII CI 
Coulomb interaction allows quasiparticles to exchange 
energy (through inelastic collisions) without coupling the 
two spin species. As shown in Figs. 0] and [5j inelastic 
scattering leaves the two distributions f+(e) and f-(e) 
strongly spin-dependent in the parallel configuration, 
while making them to coincide in the antiparallel configu- 
ration, so that both magnetoresistance and current polar- 
ization are expected to be only slightly affected. Indeed, 
TMR is only marginally affected even in the quasiequilib- 
rium regime, as shown in Fig. Illf a). where we compare 
TMR at <f> — 7T, as a function of V, for the full nonequi- 
librium and the quasiequilibrium regimes. The effect of 
energy redistribution characteristic of quasiequilibrium 
consists merely in a smoothing of some of the sharp fea- 
tures present in the nonequilibrium limit. In Fig. Illf b). 
we compare the plots of Pi at <f> — tt as functions of V for 
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FIG. 10: (color online) (a) Nonequilibrium tunnel magne- 
toresistance ratio |TMR| vs V calculated for several Y values 
at Tbath = 0.1T C , h* = 0.2A*, and <f) — tv. (b) Nonequilibrium 
Pi vs V calculated for the same values as in (a). 



both regimes. In particular, quasiequilibrium displays a 
reduction of polarization for \V\ > 2(A* — h*), and an 
increase of polarization for \V\ < 2(A* — h*). Neverthe- 
less, Pi values as large as 100% can be obtained in the 
quasiequilibrium limit as well. 

The full evolution of Pi as a function of /C co n in the rel- 
evant intervals of large polarization is shown in Fig. [TTJ 
(c-d). It is easy to notice the gradual smearing of Pj 
by increasing the electron-electron interaction strength, 
and that the quasiequilibrium limit is already reached 
for /C C oU ~ 10 (see also the discussion at the end of 
Sec. IIII Bj ). However, in FS-I-N-I-SF realistic structures 
it should be possible to keep K C al\ ~ 0.1 or smaller so 
that current polarization would be somewhat similar to 
that obtained in the full nonequilibrium limit. 

By contrast, TMR must decrease if spin- flip processes 
mix the spin-dependent distributions. In metals and at 
low temperature (typically below ~ 1 K), such processes 
are normally caused by the presence of magnetic impuri- 
ties in the N layer. Spin-flip scattering can be suppressed 
by limiting the magnetic-impurity content in the N layer, 
and by choosing <C A s f (the spin-flip relaxation length 
A s f is of the order of some /im in metals such as Cu or 
A u 49 ' 50 ). These constraints can be met fairly easily ex- 
perimentally in multilayered or planar structures like the 
ones presented here. 
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FIG. 11: (color online) (a) Tunnel magnetoresistance ratio 
|TMR| vs V in nonequilibrium (solid line) and quasiequilib- 
rium (dash-dotted line) at Tbath = 0.1T C , h* = 0.2A*, and 
cj> — tt. (b) Current polarization Pi vs V in nonequilibrium 
(solid line) and quasiequilibrium (dash-dotted line) calculated 
for the same values as in (a) . (c-d) Pi vs V in the intervals of 
large polarization calculated for several /C co ii values and the 
same parameters as in (a). 



VI. SPIN-DEPENDENT EFFECTIVE 
TEMPERATURE 

Even in the nonequilibrium case, it is interesting to 
characterize the distribution function through an "effec- 
tive" temperature and an "effective" chemical potential. 
Such effective quantities can be meaningfully defined by 
fictitiously connecting the N region to a large normal- 
metal reservoir through an insulating barrier that is suf- 
ficiently transparent to allow the flow of quasiparticlcs, 
but opaque enough not to alter the nonequilibrium con- 
dition of the electrons. One can identify the effective 
chemical potential (/i£ ff ) °f the out-of-equilibrium elec- 
tron gas in the N region with the chemical potential the 
reservoir must possess in order for the particle current to 
be zero. The effective temperature (T° B ), on the other 



hand, is taken to be equal to the one the reservoir must 
have in order for the heat current to be zero^k These 
conditions can be expressed respectively by the following 
two equations: 

J de[.U(e) -f (e- if)] = (19) 

J de e[f a (e) -f (e- fjf, Tf)] = 0, (20) 

where we have assumed that the DOS of the N reser- 
voir is equal to that of the N layer. We wish to warn 
the reader that the words "temperature" and "chemical 
potential" have to be taken in a loose sense, especially 
when the distributions are very different from equilibrium 
functions. They are merely two parameters which grasp 
important characteristic properties of the distributions, 
related, namely, to particle and heat transport. 

A general expression for the effective temperature can 
be easily derived through the Sommerfeld expansion 36 , 
obtaining 



'Cff . 



t: 



where 




dee[f <T (e)-f (e-fif,T = 0)] (21) 



de[f a (e)-6(e)\, 



(22) 



and 6(e) is the Heaviside step function. Equation Q21j) 
yields the true spin-dependent electron temperature in 
(quasi)cquilibrium. Furthermore, in the present FS-I-N- 
I-SF system, T° s turns out to depend on the strength of 
electron-electron interaction (i.e., on /C co n) as we shall 
show in the following, due to heat exchange with FS 
reservoirs with nonconstant DOS 22 ' 36 . 

In Figs. fT2Ta) and (b) we plot the calculated effective 
temperatures in the full nonequilibrium limit (i.e., for 
^Ccoii = 0), normalized with respect to the bath temper- 
ature (Tbath = 0.1T C ), as a function of the applied bias 
voltage V for = and <f> = it, respectively. In the for- 
mer case, there is no spin dependence and T|. ff is an even 
function of the bias voltage. Starting from the equilib- 
rium condition (i.e., at V = where the Tj: ff = Tbath), the 
temperature first increases reaching a maximum around 
eV/A* = 0.8 and thereafter decreases down to the min- 
imum around eV/A* — 2. The initial increase is the 
"anomalous heating" due to the presence of a finite DOS 
within the superconducting gap 22 ' 36 , while the minimum 
reflects the usual electron "cooling" which is maximized 
for voltages around twice the value of the gap22 (see the 
discussion in Sec. IIII Ajl . Such features are present also 
in the absence of an exchange field, with negligible quan- 
titative differences. By contrast, the effect of the ex- 
change field is manifest in the antiparallel configuration. 
Indeed, the effective temperatures are different for the 
two spin species, though related according to the relation 
T£ (— V) — T^(V). In particular, minima are shifted by 
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FIG. 12: (color online) (a) Nonequilibrium spin-dependent 
electron effective temperature T^ ff vs V calculated for = 
at Tbath = O.lTc and h* = 0.2A*. (b) T° ff vs V calculated for 
<j> = 7r at Tbath = O.lTc and fe* = 0.2A*. 
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an amount equal to 2ft* towards higher (lower) voltages 
for spin up (down) electrons. This fact can also be under- 
stood through the schemes shown in Fig. [7J as due to the 
shift, introduced by the exchange field, of the DOS of the 
two superconductors in opposite directions. For spin up 
electrons, maximum cooling occurs for eV — 2 (A* —ft*) 
[see Fig-Efb)] and for eV = -2(A* + ft*) [see Fig.^d)]. 
For spin down electrons maximum cooling occurs at op- 
posite voltages. The position and amplitude of maxima 
turns out to be function of the exchange field, as well as 
of the smearing parameter T. Notably, as shown in Fig. 
[T27 b). the spin-dependent effective temperatures can be 
largely different upon voltage biasing the structure. 

In Fig. ll3f a) we plot the nonequilibrium effective tem- 
perature difference ST cS = Tf - T1 B normalized to the 
bath temperature (Tbath = 0.1T C ) versus V at <fi = ir and 
for different values of the exchange field ft* . i5T eff is odd 
in the applied voltage and is more pronounced for larger 
values of ft*. For positive values of V, the maximum (as 
large as 500 % at this bath temperature) corresponds to 
the anomalous heating (occurring around eV/A* ~ 1), 
while the minimum occurs for eV/ A* > 2, and moves to 
higher values as ft* increases. 

The effect of electron-electron inelastic collisions is ad- 
dressed in Fig. fT^f b) which shows <5T off /Tb a th at cf> = ir as 
a function of the collision strength IC co \i for eV = 1.0 A*, 
ft* = 0.2A*, and T ba th = 0.1T C . We find that a dra- 
matic effect of electron-electron interaction, that leads 
to a strong suppression of the effective temperature dif- 
ference on the scale of /C co ii — 0.1. With the material 
parameters given in Sec. IIIIBI this would correspond to 
a N region with « 470 nm. 

The possibility to have access to different spin- 
dependent electronic temperatures suggests that we in- 
vestigate the potential for the implementation of spin- 



FIG. 13: (color online) (a) Difference of the nonequilibrium 
spin-dependent electron effective temperatures ST eS vs V cal- 
culated for different h* values at Tbath = 0.1T C and (j> — it. 
(b) 5T cB vs JCcM calculated at eV = 1.0A* for T ba th = 0.1T C , 
h* = 0.2A*, and <j> = n. 

tronic devices. In particular, we conceive a nanostruc- 
ture like that shown in Fig. I14f a) where an additional 
voltage-biased superconducting electrode (S3) is tunnel 
coupled to the N region through a junction of resistance 
R p » Rir^, while the exchange fields in Si, 2 are arranged 
in the antiparallel configuration (cf> = ir). The presence 
of a superconducting extraction lead is crucial, since in 
the tunneling process the quasiparticle current through 
S3 will depend exponentially on the electron tempera- 
ture in On one hand, the setup considered al- 
lows direct measurement of the spin-dependent electron 
temperatures and S3 may act as a thermometer—. S3 
provides, in fact, access to the whole distribution func- 
tions f a (e) from the voltage-dependent differential con- 
ductance of the NIS3 junction^. On the other hand, upon 
biasing the S3 electrode with T4rias, the existence of dif- 
ferent spin-dependent temperatures in the normal metal 
region yields a finite current polarization Pj defined in 
the usual way as 

p ,„ v , _ I S + 3 (V,V^s,h*)-I s J(V,V hias ,h*) 

I{ ' Uas ' } - il>(v,v hias ,h*) + i s j(v,v hias ,h*y 

(23) 

where 

4 S3 (^bia S ,ft*) =4" / d£ ^ S3 ( £ ) 

eR p J (24) 
x[fa(e,V,h*)- f (e + eV hias )}, 

and J\f Ss (e) is the normalized DOS of S3 . In the following 
we assume for simplicity that Af 3 (e) is identical to the 
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FIG. 14: (color) (a) Sketch of a possible setup for the imple- 
mentation of a spin-polarized current source. An additional 
superconducting electrode (S3) is coupled to the N region 
through a tunnel junction of resistance R p . The exchange 
fields in Si,2 are arranged in the antiparallel configuration. 
Spin-polarized current can be extracted by biasing the S3 ter- 
minal with 14>ias- (b) Contour plot of the nonequilibrium cur- 
rent polarization Pi vs V and Vbias at Tbath = 0.1T C and 
h* = 0.2A*. (c) Pi vs JC coU calculated eV = 1.0A* and 
e Vbias = 0.5 A* for T bath = 0.1T C , h* = 0.2 A*, and <b = n. 



density of states of Si^ in the absence of an exchange 
field (i.e., h*=0). 

Figure fT4Tb) shows the calculated nonequilibrium Pj 
as a function of V and Vbias, for Tbath = 0.1T C and 
h* = 0.2 A*. For easily attainable values of V and Vbias 
pure (i.e., 100%) spin-polarized current of both the par- 
allel and antiparallel spin species can be achieved. Fur- 
thermore, we note that |P/| largely exceeds 50% over a 
wide region in the (V, Vbias) plane. It is worthwhile to 
note that replacing the S3 terminal with a N metal reser- 
voir would completely hinder the extraction of a finite 
spin-polarized current for <f> = it. This stems from the 
insensitivity of the tunneling current in a NIN junction 
to the N region temperature. The <f> = case for a N lead 
connected to a similar setup was analyzed in Ref. [3(| 

The role of electron-electron inelastic relaxation is dis- 
played in Fig. ITHc) which shows Pi versus K ca \\ at 
eV = A* and eV bim = 0.5 A*, for Tbath = 0.1T C and 
h* = 0.2 A*. We find that the suppression of the cur- 



rent polarization occurs on a larger range of values of 
^Ccoii with respect to electron effective temperatures dif- 
ference [note the different scale for the horizontal axis 
of Fig. rnif b) and fT4T c)]. At K, co \\ = 0.1, for example, 
Pi is reduced only by about 35%. This behavior is not 
surprising: it originates, as mentioned above, from the 
fact that the current extracted from S3 is exponentially 
sensitive to T° ff . As a consequence, even a small tem- 
perature difference between spin species yields a large 
current difference, and gives rise to a sizable Pj vs /C co n 
characteristic. 

We shall further comment the impact of relaxation in 
this setup as compared to that of Fig. 1 by inspecting 
Fig. 14(c) and Figs, ll(b-d). In particular they show 
that while electron-electron interaction is crucial in sup- 
pressing Pi in the former case, it is much less important 
in the setup of Fig. 1. The reason stems from the fact 
that while in the present case the electric current is spin- 
dependent only thanks to presence of spin-dependent dis- 
tribution functions in the N region (as expressed by Eq. 
24), in the setup of Fig. 1 the spin selectivity originates 
from both the distribution functions f a and the spin- 
dependent superconducting DOS (see Eqs. 10 and 11). 
From this follows that, in the first setup, current polar- 
ization will completely vanish at quasiequilibrium (where 
the quasiparticle distribution functions result to be iden- 
tical and spin-independent), while in the second system 
spin polarization will persists also for identical thermal 
distribution functions owing to the additional spin selec- 
tivity provided by spin-split DOS in the superconductors. 



VII. CONCLUSIONS 

In this paper we have analyzed the nonequilibrium 
spin-dependent transport properties in superconductor- 
normal metal tunnel nanostructures, where the supercon- 
ductors present a proximity-induced effective exchange 
field (equivalent to that of a superconductor in a mag- 
netic field). The latter is due to the proximity of a fer- 
romagnetic thin layer which exerts, under appropriate 
conditions, a nonlocal influence. We have computed the 
quasiparticle distributions of a normal metal layer tunnel 
coupled to two superconductors with non-collinear ex- 
change fields. In the full nonequilibrium limit, i.e., in the 
absence of any inelastic relaxation mechanisms, the dis- 
tribution functions depend on spin and display unusual 
features, such as population inversion and double-step 
shape, depending on the bias voltage applied to the su- 
perconductors. Spin-dependence persists also in the pres- 
ence of inelastic Coulomb interaction, which produces a 
smoothing of the sharp features and tends to drive elec- 
trons into the quasiequilibrium regime, where the energy 
relaxation is strong enough to allow electrons to thermal- 
ize. This interesting behavior is reflected in the current- 
voltage characteristic, which shows a strong dependence 
on the relative orientation of the exchange fields existing 
in the superconductors. Notably, even in the quasiequi- 
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librium regime, a huge tunnel magnetoresistance ratio 
and a complete current spin-polarization were found over 
a wide range of bias voltages and for realistic parameters. 
The impact of the exchange field as well as of the pres- 
ence of subgap states in the DOS of the superconductors 
have been addressed. 

We have characterized the out-of-equilibrium distribu- 
tion functions through an "effective temperature" and 
an "effective chemical potential" , defined in an opera- 
tive sense. We have found that such effective temper- 
atures are strongly spin-dependent for anti-parallel ex- 
change fields, the relative temperature difference being 
as high as 500%. Furthermore, we have discussed the 
possibility of producing spin-polarized currents by cou- 
pling the N region to an additional superconducting lead, 
finding that 100% spin-polarization is realistically achiev- 
able. This effect is fairly robust against the occurrence 
of inelastic electron-electron collisions. 

We shall finally comment on some possible applica- 
tions of the structures here presented. An immediate 
first application of this system is the implementation 
of storage cell elements, thanks to the very large TMR 
values [see Fig. Ufa)]. Magnetic-field-controlled current 
switches can be envisioned as well [see Fig. [Ha)]. Im- 
portantly, power dissipation is intrinsically limited ow- 
ing to the small currents driven through NIS junctions. 
For example, assuming R t — 10 3 f2 and aluminum (Al) 



electrodes at T — 0.1T C « 0.12 K, a dissipated power 
in the range of 10~ 15 . . . 10~ 12 W can be achieved for 
2(A* - h*)/e < \V\ < 2A*/e. This makes this setup 
attractive for low-dissipation cryogenic applications. In 
light of a realistic implementation, ferromagnetic alloys 
such as Cui-xNia^ 3 - or Pdi-zNL^ 4 - (which allow fine tun- 
ing of h through a proper choice of x) are promising can- 
didates. For example, in Pdi-^Ni^ alloy with x = 0.1, 
h ~ 10 meV resulting in £p ~ 5 nm 54 . By choosing 
Al as S electrodes (with A ~ 200 [ieV and £s » 300 
nm£) it turns out that h* in the range ~ 0.2A*...0.5A* 
can be achieved. In such nanostructures the bias voltage 
con be fed through outer normal metal electrodes, tunnel 
coupled to the ferromagnetic layers in order to prevent 
depolarization of the F electrodes. This would result in 
adding in series extra resistances which could be easily 
engineered in order to minimize the induced correction 
to both tunnel magnetoresistance and operating voltage. 
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